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Abstract 

"Acoustic spacetimes" , in which techniques of differential geometry are 
used to investigate sound propagation in moving fluids, have attracted consid- 
erable attention over the last few decades. Most of the models currently con- 
sidered in the literature are based on non-relativistic barotropic irrotational 
fluids, defined in a flat Newtonian background. The extension, first to special 
relativistic barotropic fluid flow, and then to general relativistic barotropic 
fluid flow in an arbitrary background, is less straightforward than it might at 
first appear. In this article we provide a pedagogical and simple derivation of 
the general relativistic "acoustic spacetime" in an arbitrary (d+1) dimensional 
curved-space background. 
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1 Introduction 

In this article we shall present a simple and pedagogical derivation of the general 
relativistic version of the "acoustic metric" defined on an arbitrary curved (d + 1) 
dimensional background spacetime. While there are related observations and more 
limited derivations extant in the literature, we feel that the current analysis has 
some definite advantages. For instance, the early 1980 analysis due to Moncrief is 
restricted to perturbations of a spherically symmetric fluid flow on a Schwarzschild 
background [I], and in the more recent 1999 derivation due to Bilic [2] it can be 
somewhat difficult to discern what is truly fundamental input from what is derived 
output. Since scientific interest in this field is both significant and ongoing [3j SJ [5], 
we feel it useful to carefully lay out the minimal set of assumptions and logic flow 
behind the derivation. 
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Our strategy is as follows: 

• We shall first motivate the result (up to a conformal factor) by considering the 
acoustic version of the Gordon metric, which was introduced for geometrical 
optics in 1923 [6]. 

• We shall carefully specify what we mean by barotropic flow, the assumption 
that energy density is a function of pressure only. (Specifically, barotropic flow 
is more general than either isothermal or isentropic flow.) 

• We shall then carefully specify what is meant by "irrotational flow" in a general 
relativistic context, introducing the appropriate notion of velocity potential. 

• From the relativistic Euler equation, using only the irrotational condition and 
the barotropic condition, we will derive the relativistic Bernoulli equation. 

• From the relativistic energy equation, using only the barotropic condition, (that 
is, without using the irrotational condition), we shall derive a flux conservation 
law (continuity equation). 

• We shall delay the introduction of thermodynamic arguments as long as practi- 
cable. (We would argue that thermodynamics is in fact a side issue not central 
to the derivation.) 

• As usual, the acoustic metric follows from combining the linearized Bernoulli 
equation and linearized equation of state with the linearized continuity equa- 



• We shall carefully explain the subtleties involved in taking the non-relativistic 
limit. 

• We shall finish with some discussion, and relegate several thermodynamic ob- 
servations to the appendices. 

The key result that we shall be aiming for is this: the (contravariant) acoustic metric 
governing acoustic perturbations of an irrotational barotropic fluid flow in (d + 1) 
dimensions is 



tion. 




(1) 
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In counterpoint, the (covariant) acoustic metric in (d + 1) dimensions is 



Gab 




(2) 



Here: 

• c s is the speed of sound, defined as usual via cj. = c 2 dp/dg. 

• c is the speed of light, used for instance in defining x° = ct. 

• V is the dimensionless 4- velocity of the background fluid flow. 

• h ab = g ab + [Vo] a [Vo]& is the dimensionless orthogonal projection of the physical 
spacetime metric g a b onto the 3-space perpendicular to the 4-velocity of the 



• the indices on the background fluid flow [V ] a are lowered and raised using the 
physical spacetime metric g ab and its inverse g ab . 

• in contrast Gab and G ab are defined to be matrix inverses of each other; these 
indices are not to be raised and lowered with the physical spacetime metric. 

• no is the background number density of fluid particles. 

• £>o is the background energy density of the fluid. 

• po is the background pressure of the fluid. 

If one appeals to relativistic thermodynamics, by introducing the specific enthalpy 
w = (g + p)/n, then, as discussed in the appendix, the acoustic metric can be 
somewhat simplified: 



These key results are easy to motivate (but not derive) in the limit of "ray acoustics", 
also known as "geometric acoustics" , where we can safely ignore the wave properties 



fluid. 




(3) 
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of sound. In this limit we are interested only in the "sound cones" . Let us pick a 
point in spacetime where the background fluid 4- velocity is Vq. Now adopt Gaussian 
normal coordinates, and go to the local rest frame of the fluid. Then taking x° = ct 
we have 

[V ]°->(i; 0,0,0), 



and 



9ab -»■ 



-10 
10 
10 
1 



h a b — dab 



a [V0\b ->• 





10 

10 

1 



In the rest frame of the fluid the sound cones are (locally) given by 

-<?„ dt 2 + \\dx\\ 2 = 0, 



(5) 



(6) 



(7) 



which we can rewrite as 

2 2 

- % (cdt) 2 + ||d£|| 2 = -% (dx ) 2 + ||df|| 2 = 0, 

implying in these special coordinates the existence of an acoustic metric 

-c 2 /c 2 

10 

10 

1 



Gab OC 



That is, transforming back to arbitrary coordinates: 



(8) 



(9) 



G ab oc — § [V c 



a \ V0 b 



+ h 



ab ■ 



(10) 



We now rewrite this as 



Gab^-^[V }a[V } b + {g ab 



Oja 



o}b} oc g ab 



c 2 

1 - ^ 

c 2 



a \VQ\b 



(11) 



Note that this is essentially a generalization of the derivation (dating back to 1923) of 
the so-called "Gordon metric" [6] used in ray optics to describe the "optical metric" 
appropriate to a (possibly) relativistic fluid with position-dependent refractive index 
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n(t,x). In Gaussian normal coordinates co moving with the fluid the Gordon metric 
is given by 

-1/n 2 
10 
10 
1 



Gab oc 



(12) 



That is, transforming back to arbitrary coordinates: 



Gab OC 



o] a [V )b + {g a b + [V )a [V }b} oc g ab 



a [VO b 



(13) 



Note that in either the ray acoustics or ray optics limits, because one only has the 
sound cones or light cones to work with, one can neither derive nor is it even mean- 
ingful to specify the overall conformal factor [31 El E]. (See also reference [10] for an 
equivalent Godon-like metric based on relativistic fluid dynamics but without calcu- 
lation of the overall conformal factor.) The calculation presented below is designed 
to go beyond the ray acoustics limit, to obtain a relativistic wave equation suitable 
for describing physical acoustics — all the "fuss" is simply over how to determine the 
overall conformal factor (and to verify that one truly does obtain a d'Alembertian 
equation using the conformally fixed acoustic metric). 



2 Basic general relativistic fluid mechanics 

In this Section we consider fully relativistic barotropic inviscid irrotational flow on an 
arbitrary general relativistic background, and derive the relevant wave equation for 
linearized fluctuations. We already know exactly what happens in the non-relativistic 
case [31 [7J EJ [9] . We go straight to curved spacetime relativistic fluid mechanics, where 
the geometry is described by a metric tensor g a b(x) of signature — h ++, [or more 
generally — (+) d ], and the fluid is described by the energy density g, pressure p, and 
4-velocity V a where 

g a bV a V b = -l. (14) 

The two relevant fluid dynamical equations are extremely well-known (see, for in- 
stance, Hawking and Ellis [11J). 

Relativistic energy equation: 

V a (gV a )+p(V a V a ) = 0. (15) 
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Relativistic Euler equation: 

(f? + p)V h V b V a = -(g ab + V a V b ) V b p. (16) 

These equations can be combined into the single statement that the stress-energy 
tensor is covariantly conserved (see, for instance, Hawking and Ellis [TT] ) 

V a [(Q + p)V a V b +pg ab ] = 0. (17) 

Note the subtle (perhaps not so subtle) differences from the non-relativistic case: g 
is now energy density, not mass density p, and the continuity equation looks different 
- at least at this stage of the calculation it has been replaced by an "energy con- 
servation" equation, and the Euler equation now contains the combination (g + p). 
The 4-acceleration of the fluid is 

A a = V b V b V a , such that A a V a = 0. (18) 

2.1 Defining relativistic barotropic flow 

The fluid is barotropic if the energy density is a function only of the pressure - 
that is, there is some function such that g = g(p). Appealing to the inverse function 
theorem we will have p = p(g) at least piecewise on suitable open sets in energy 
density. This barotropic condition is much more general situation than assuming 
isothermal or isentropic flow. Specifically, for a one-component fluid: 

• Any zero-temperature fluid is automatically barotropic. 

• Any non-zero temperature but isothermal fluid is automatically barotropic. 

• Any zero-entropy fluid (superfluid) is automatically barotropic. 

• Any isentropic fluid is automatically barotropic. 

However, merely assuming the barotropic condition does not imply any of the isen- 
tropic, zero-entropy, isothermal, or zero-temperature conditions. 

2.2 Defining relativistic irrotational flow 

What do we mean by irrotational flow in a relativistic setting? Construct the 1-form 

V = V a dx a , (19) 
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and consider the 2-form 002 and 3- form t^. 



U2 = dv; t 3 = v A U2- (20) 

In a relativistic context, setting the 2-form 0J2 = is too strong a condition, setting 
the 3-form t 3 = is just right. (See, for instance, Hawking and Ellis [H].) Indeed, 
adopting Gaussian normal coordinates and going to the local rest frame of the fluid, 
where V a — > (1; 0, 0, 0), setting £3 = implies duVj] = 0, so the spatial components of 
the flow velocity are locally irrotational (in the sense of being curl-free). But £3 = 
implies (via the Frobenius theorem) that locally there exist functions a and such 
that 

v = adQ; V a = a g ab V fe 0. (21) 
But then, from the normalization condition for the 4-velocity 

V = - g " V * 9 (22) 

V-9°" v„e v„e 

We shall find it extremely useful to define 

||V0|| 2 = -g ab v a e V b 6 (23) 

so that 

V = Qf. (24) 

This is the relativistic condition for irrotational flow. The function can now be 
interpreted as the general relativistic version of the velocity potential. Note that for 
any smooth function F(-) we can replace -H- F(Q) without affecting the 4-velocity 
V a — both numerator and denominator above pick up factors of F'(Q) which then 
cancel. This freedom in choosing the scalar potential, 0, will be very useful when 
analyzing the Euler equation, and will allow us to obtain a relativistic Bernoulli 
equation. 

For the discussion below define the projection operator 

h ab = g ab + ya yb _ ^ 

Using this projection operator, and the 4-orthogonality of 4-velocity with 4-acceleration, 
it is easy (for relativistic irrotational flow) to calculate 

V b V a = (5 a c + V a V c )^-. (26) 
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Therefore, for relativistic irrotational flow, the 4-acceleration reduces to 

A a -VV b V a - ^6 a +V a V) _^ Va@Vb@ , (27) 

so that 

A a = -(^ + V ro V r6 )V 6 (log||Ve||). (28) 

Note that automatically A a V a = 0, as required. Furthermore, note that the 4- 
acceleration is the projection of the gradient of a scalar, and that can be chosen 
to have the dimensions of a distance so that VO is dimensionless. 



2.3 From Euler equation to Bernoulli equation 

The Euler equation for relativistic irrotational flow now reads 

(g ab + V a V b ) V 6 (log 1 1 Ve||) = (g ab + V a V b ) (29) 

Q + P 

We now make use of the barotropic condition g = g(p) to obtain 

VbP - fP dp (30) 



Q + P Jo Q(p)+P' 
so the Euler equation becomes 

(g - + v^)v t (->og||ven + j ( P ^)=o. pi) 

Note that for any arbitrary function /(©) we have 

(^ + V -y 6 )VJ(0)=O, (32) 

(since the projection operator kills the gradient). Therefore (using both the irrota- 
tional and barotropic conditions) we can integrate the Euler equation to yield: 

f p dp 

— log ||Ve|| + / — + /(6) = constant, (33) 

Jo QVP) + P 

where /(©) is (at this stage) an arbitrary "function of integration". This is our 
preliminary version of the general relativistic Bernoulli equation. To see how we 



9 



might further simplify this result, recall (see for example [U IU [9]) that for non- 
relativistic irrotational inviscid barotropic flow the non-relativistic Bernoulli equation 
is 

1 f p do 

6 + -(V0) 2 + / -A + f(t) = constant. (34) 

2 Jo Pip) 

(In the non-relativistic case p is the mass density.) Now, in the non-relativistic case 
we can always redefine 

0^0+/ f{t)dt, (35) 



and use this transformation to eliminate f(t) — so we can without loss of generality 
write 

, r dp 



+ -(V0) 2 +/ -r-r = constant. (36) 
2 Jo Q{P) 

In the relativistic case we now note: 

-log 1 1 ve|| + /(e) = -io g ( e - /(e) ||ve||) = -io g ||VF(e)||. (37) 

That is, in the relativistic case, by making the transformation 

r e 

F(S) = / e" /(e) d9, (38) 



we can (without changing V a ) absorb the arbitrary function /(■) into a redefinition 
of 6, and so write 

r p dp 

- log 1 1 V6|| + / — -v- — = constant. (39) 
Jo Q(P) + P 

Finally, without loss of generality we can rescale Q to set the constant appearing 
above to zero and so obtain: 

Jog II ve|| = [ P -^—- (40) 

Jo Q{P) + P 

This is our final form for the general relativistic Bernoulli equation. Note that 
there is no longer any freedom left in choosing 0, we have now used it all up. It 
is relatively common to exponentiate the above and rewrite the general relativistic 
Bernoulli equation as 

||V9|| =exp ( [ P . dp V (41) 



'o Q(p)+P, 

For the sake of pedagogical development of the argument we have not yet rewritten 
the integral on the RHS in terms of other thermodymamic variables; we prefer to 
delay this for now. 
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2.4 From energy equation to flux conservation 

The energy conservation equation is 

v a (f?n+p(v a n = o, 

which becomes 

V a V a g+(g + p)(V a V a )=0, 



or 



Q + P 



(42) 

(43) 
(44) 



Using the barotropic condition g = g(p), (but note, now without using the irrotational 
condition), this can be written as 



V a V a 



j« dg 

l e(p=0) q + p(q) 



+ {V a V a ) = 0, 



(45) 



which implies 

V a V a \ exp 



dg 



8( P =0) 



Q + p{o) 



cxp 



dg 



0(p=O) 



(V a V a ) = 0. (46) 



This can now be rewritten as a continuity equation (flux conservation equation): 

dg 



V« < exp 



0(p=O) 



q + p{q) 



ya 



0. 



(47) 



This is now a "standard form" zero-divergence continuity equation — note this the 
ability to derive this transformed form of the energy equation depends only on the 
barotropic assumption. Also note that we have not yet needed to even introduce, 
let alone discuss in any detail, the particle number density n, nor introduce any 
thermodynamic arguments. 



2.5 Introducing the number density 

Subject only to the barotropic condition, we have just derived the flux conservation 
equation given immediately above. We now suppose the barotropic fluid contains 
some type of conserved "tracker" particles. For example, one might be interested 
in counting baryons, Fe 56 nuclei, or leptons. We now explicitly assume translation 
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invariant and time invariant composition of the fluid. (That is, we assert that the 
equation of state is the same throughout the fluid; in other words there is no ex- 
plicit time or position dependence in the equation of state.) Therefore the ratios 
of these tracker particle densities must be translation and time invariant constants. 
Furthermore since these "tracker particles" are all assumed to be conserved 



V a {n 4 V a } = 0, 



(46 



while, since the fluid is assumed to be barotropic, there must be functions of pressure 
p and energy density g such that 



rii = m(p) = rii(g). 
But the only way to satisfy all these constraints is if 



Ui(p) 



n 



i(p=0) ex P 



> Q(p=0) Q + p(g) 
In particular for the total particle density we have 

•<?(p) d £ 



n(p) 



n 



(p=o) exp 



Q(p=o) 



+ p(g) 



(49) 



(50) 



(51) 



This observation is extremely convenient in that allows us to physically interpret the 
quantity 



cxp 



d^> 



e(p) 

e(p=o) Q + P{0) 



n(p) rii(p) 



(52) 



n(p=o) nj(p=o) 

as being proportional to the number density of constituents making up the fluid. (If 
one wishes to take an extreme point of view and eschew all thermodyamic arguments 
completely, one could simply take this equation as the definition of a "shorthand 
symbol" n(p), and ignore the physical interpretation of n(p) as particle number 
density.) In terms of the number density the conservation equation now simply 
reads 

V a {nV a } = 0. (53) 

Furthermore, note that using these results we can rewrite the relativistic Bernoulli 
equation (T4T|) as 

rp d[g(p)+p] f 8 dg 



log ||ve| 



dp 
g(p) +p 



g(p) + p 



log 


~g(p) +p~ 


- log 


n(p) 


= log 




Q(p=o) 









e (P =o) Q + p(q) 
[g(p) + p] n (p=0) 
nip) Q(p=o) 



(54) 
(55) 
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That is 



3 Linearization 



Let us now write 



ivei 



[g(p) + p] n ( p =0) 
n{p) Q(p=o) 



(56) 



e = e + eei + 



which in particular implies that 



and further assert 



V = V + eV 1 + . 

g = Q + e g 1 + . 
P = Po + epi + • 



(57) 

(58) 

(59) 
(60) 



Using these relations we now linearize the fluid equations around some assumed 
background flow. (Note that both the background fluid flow (Vo, f?o, Po), and the 
linearized fluctuations, satisfy the Bernoulli and energy conservation equations.) In 
a wider context, extending far beyond fluid dynamics, we mention that it is quite 
common for linearized fluctuations around an appropriately defined background to 
exhibit an "effective spacetime" behaviour [31 [121 US] ■ 

3.1 Linearized continuity equation 

From the continuity equation we see 

"0O+££>1 + -- 



V a jexp 

Then to first order in e 
V a <^ exp 



dg 



e(j,=o) 



on 



q + p(q) 



K + eV? + 



(61) 



Qi 



dg 

c(p=o) 6 + p{q)] \6o + P« 
Using the number density we can rewrite this as 

Ql 



V a + 



(62) 



V a < n(p ) 



V Q a + 



Qo+Po 



0. 



(63) 
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3.2 Linearized irrotational flow 

Perturbing the 4-velocity for relativistic irrotational flow yields 

g ab V 6 (0 O + e6i + . . . ) 



V + eV 1 + 



y/-g°*> V a (0 o + e0! + . . . ) V 6 (0 O + eB 1 + ...)' 

Then expanding to first order in e we see 

(g ab + V a V b ) V b 0! 
1 " l|V0o|| 

Now use the relativistic Bernoulli equation to write 



(rpo 



dp 



Alternatively 
Note that we automatically have 



q(p) + V 



v? = . n " 0( ' =o) {g ab + V«V b ) v 6 e!. 

[00 + PoJ ^(p=0) 



^K a V^ = 0, 



as required by the normalization condition, g a bV a V = —1, for V a . 

3.3 Linearized equation of state 

Linearizing the equation of state we see 

Po + epi H = p(Qo + eg ± + ...), 

that is 

dp 

Qi, 



dp 
dg 



80 

which we use to define what we shall soon enough see is the speed of sound 

pi = -o ei- 

c 2 
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3.4 Linearized Euler (Bernoulli) equation 

Linearizing the Bernoulli equation requires us to consider 

•(po+epi+... ) 

o loo r -a''"' V ..((-)„ f rB, + ...)V,[f)„irH,l...)l / 
Then to first order in e 



1 r(po+e Pl + ... ) j 

- log [-<7 cd V c (6o + e0i + . . . ) V d (0 o + eGx + ...)] = / "TTj - ' (72) 

-^VceiVdGo pi 



-<7 cd v c e v d e ^o + po' 

That is, using the linearized equation of state, 

l|ve || c 2 ^o + po' 



(73) 



(74) 



which we can rearrange to 



c 2 s ||V0o| 



0i = -(0o + Po)^^^T- (75) 



That is, now using the general relativistic Bernoulli equation, 

Qi = -(f?o +Po)^ exp (- [ P ° ^ P ) ^ a V a 6i. (76) 
Finally, using the number density we can rewrite this as 

Qi = -~r^-- 2 Q(p=o) V" V o e!. (77) 

3.5 Deriving the d'Alembertian equation 

Now combine these results: insert the linearized Euler (Bernoulli) equation ( 177ft . 
and the linearized irrotational condition (167p . into the linearized continuity equation 
(1631). We obtain 

nl ° 2 eip=0) -VfVfr&i + n °f (p=0) A g ab + V a V b ) 1 = 0. 



n(p=o) c 2 s (Oo + Po) n(p= ) (00 + Po) 

(78) 

But of course £( p= o) and ri( p=0 ) are irrelevant position-independent constants, so we 
can just as easily write 



22 2 



VfV'VbB! + -, ^(g ab + W) V 6 8i = 0. (79) 



c 2 s (Qo+Po) (Qo+Po) 
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Introducing the projection tensor h ab = g ab + V Q a V b , and factorizing, this becomes 



06 



lie 



Qo+Po 



-3 V a V b + h 



,ab 



0. 



(80) 



This is in fact exactly the result we want — a d'Alembertian equation for the per- 
turbation in the velocity potential G. 

• After some trivial notational changes, this agrees (where the formalisms over- 
lap), with both the observations of Moncrief p], and with the presentation of 
Bilic pj, but the present exposition gives much more attention to the underly- 
ing details and makes only a bare minimum of technical assumptions. 

• Note that everything so far is really dimension- independent, and that we can 
now read off the acoustic metric simply by setting 



-QQ 



ab 



— V a V b 
-2 v o v o 



h 



ab 



Qo+Po 

The dimension-dependence now comes from solving this equation for Q ab 



(81) 



4 The general relativistic acoustic metric 

From the dimension-independent result above we have, in (d+1) dimensions, 



(_l)(d+l)/2 g(d+l)/2-l _ _ C _ 



n? 



c; \Qo + Po 



(d+1) 



(82) 



whence 



_2 \ (d+l) 



£ ( *g ) 



(83) 

ct \Qo+Po,' 

Dropping an irrelevant overall constant factor of c -2 '^ -1 ), we finally have the (con- 
travariant) acoustic metric 



tab 



„2 r -l \ -2/(d-l) 
u c s 



.Qo+Po 

and (covariant) acoustic metric 



-^V a V b + h ab 



M) 



fab 



2 r -l \ 2/(d-l) 



Qo+Po 



--§ [VoUVoh + hab 



(85) 
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If one wishes to go to the extra trouble of making the acoustic metric dimensionless 
it is easy to re-insert appropriate position-independent constants in the conformal 
factor and obtain 

/ 2 \-2/(d-l) 2 

G ab = 2 n °^ C - C - V X + h*\, (86) 

and 

/ 2 \2/(d-l) 2 

Gat = 2 n ° g(p =° )C -% [V ] a [V ] b + h ab \. (87) 

Of course these extra position-independent constants in the conformal factor carry 
no useful information and are commonly suppressed. 



5 The non-relativistic limit 

Compare this general relativistic acoustic metric with the non-relativistic limit, 
where the coordinates are most conveniently taken to be x a = (t; x l ) and where 
the d'Alembertian equation reduces to [3j EJ |9] 



. V a V b -\- h ab 

c 2 V V + 11 



0. 



with po 7^ now being the non-relativistic mass density, (not the relativistic energy 
density), and the meanings of Vq and h are suitably adjusted. In the non-relativistic 

case 



Kf = (l;n h ab 
and independent of dimensionality we have 

V^G G ab = po 

implying 



" 








6 ij 



1 



— —V a V b -I- h ab 
2 V V + a 



->ab 



which can be inverted to yield 



lab 



Po\ 

C s J 



-2/(d-l) 



2/(d-l) 



;V a V b + h ab 



1) 


\-(c 2 s -V 2 ) 


-v j 






5 ij 



19) 



(90) 



(91) 



(92) 
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To see how and under what situations the general relativistic acoustic metric 
reduces to this non-relativistic acoustic metric first consider the conformal factor: 
In the non-relativistic limit po <C go and go ~ fh no, where rh is the average mass 
of the particles making up to fluid (which by the barotropic assumption is a time- 
independent and position-independent constant). So in the non-relativistic limit we 
recover the standard result for the conformal factor [3j El [8j |9] 



"6 c s 



no 



1 Po Po 
oc — . 



m 2 c 



(93) 



go + Po mc s 

To now probe the tensor structure of the non-relativistic limit, let us recall that: 

• c is the speed of light. 

• c s is the speed of sound. 

• v is the three-velocity of the fluid. 

Take conventions so that the physical spacetime metric and four-velocity are both 
dimensionless. In particular, the coordinates are chosen to be x a = (ct; x l ). Now 
adopting Gaussian normal coordinates at the point of interest 



g a b = Vab = 
For the four-velocity we in general have 

^ = 7(i;«*/c); 



-10 
10 
10 
1 



(94) 



Ob 



7(-l;t//c) 



(95) 



Remember that dx° = c dt, and note that the 7 factor is defined using c the physical 
speed of light. 

Now let us take the nonrelativistic limit. Ignoring the conformal factor, which 
comes along for the ride, we have 



Gab oc 



9ab + 1 - 



oWoli 



Then for the time-time component of the acoustic metric 

,,2 - 



£00 oc 



■1 + 



1-^ 

c 2 



7 



(96) 



(97) 
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This implies 



Goo oc 



c 2 s -v 2 , Q(y\cy) 



c 4 



(98) 



So as long as both the speed of sound and the 3-velocity of the fluid are small 
compared to the speed of light we are justified in approximating 



c 2 -v 2 
Goo oc — 



+ 



(99) 



In contrast for the time-space components of the acoustic metric 



Goi oc 



1-3 7 2 ("1) 



This implies 



£(H oc 



0(clv 2 )\ v 



(100) 



(101) 



So as long as 6ot/i the speed of sound and the 3-velocity of the fluid are small 
compared to the speed of light we are justified in approximating 



Goi oc h 

c 



(102) 



• Finally for the space-space components we note 



Gij oc 



c: 



This implies 



Gij oc [c% + 0(v 2 /c 2 )] . 



(103) 
(104) 



So as long as the 3-velocity of the fluid is small compared to the speed of light 
we are justified in approximating 



Gij oc 5ij + 

Collecting these results we see that in the nonrelativistic limit 



(105) 



Gab OC 



■-(c 2 -v 2 )/c 2 


—v l jc 




Sij 



+ ... 



(106) 
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But now we realise that in the present context c is just some convenient fixed con- 
version constant from dx° = cdt to dt, so if we work in terms of the coordinates 
(t; x l ), which are perhaps more natural in the non-relativistic limit, then 



oc 







-V 1 


5ij 



+ 



(107) 



as required. Note that to do all this it is essential that both v and c s are small 
compared to c, though there is no constraint on the relative sizes of v and c s . 



6 Discussion 

Under what conditions are the fully general relativistic derivation of this article 
necessary? (The non-relativistic analysis of [31 [TJ El [9] is after all the basis of the bulk 
of the work in "analogue spacetimes", and is perfectly adequate for many purposes.) 
The current analysis will be needed in four separate situations: 

• When working in an arbitrary curved relativistic background; 

(for example in the problems considered by Moncrief PQ, and Bilic [2]). 

• Whenever the fluid is flowing at relativistic speeds; 

(for example in the problems considered by Moncrief PQ, and Bilic [2]). 

• Less obviously, whenever the speed of sound is relativistic, even if background 
flows are non-relativistic; 

(for example a near-equilibrum photon gas where c 2 s = |c 2 but flow velocities 
are all small v <C c). 

• Even less obviously, when the internal degrees of freedom of the fluid are rel- 
ativistic, even if the overall fluid flow and speed of sound are non-relativistic. 
(That is, in situations where it is necessary to distinguish the energy density g 
from the mass density p; this typically happens in situations where the fluid is 
strongly self coupled — for example in neutron star cores [14] or in relativistic 
BECs @].) 

In developing the current derivation, we have tried hard to be clear, explicit, and 
minimal — we have introduced only the absolute minimum of formalism that is re- 
quires to do the job, and have eschewed unnecessary side issues. We hope that the 
formalism will be useful to practitioners in the field of "analogue spacetimes" , partic- 
ularly with regard to ongoing and future developments [4J |5]. In particular, even in 
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the non-relativistic case it is already known that adding vorticity greatly complicates 
the situation [15], and a deeper general relativistic analysis of this situation would be 
interesting. Looking further to the future, the "fluid-gravity correspondence" hints at 
even deeper connections between curved spacetimes and fluid dynamics (161 [T71 [18]. 
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Appendix: Specific enthalpy and thermodynamic 
considerations 

We now consider some thermodynamics which for pedagogical purposes we have 
delayed as much as possible. (Though delayed, we emphasize that thermodynamics 
does have important consequences — see for instance [15]-) Suppose we take the 
specific enthalpy as primary 

g + p 



w 



n 



and use the fact that we have already deduced 



n = n (p= o) exp 

to then write 

w = exp[log(£> + p)]n~ = exp 
This implies 



dg 



Q( P =o) 



d[g + p\ 
g + p 



exp 



dg 



Q(p=0) 



q + p(q) 



Q(p=o) 
w = — — - exp 

^(p=0) 



p dp 



o Q + P. 



( P =o) exp 



dp 
o Q + P. 



(108) 



(109) 



(110) 



(HI) 



This so far is a purely (barotropic) thermodynamic result, independent of any irro- 
tational condition. 
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A secondary result, now specifically tied to the Bernoulli equation (and hence to 
irrotational flow), is that 

w = W(p=0} ||V9||. (112) 

Note that the way we have presented the derivation we have been able to delay and 
avoid the need for thermodynamic arguments as far as possible. (In contrast, this 
equation is the starting point adopted by Bilic in his derivation of relativistic acoustic 
geometry [2], what for us is a peripheral result has in that analysis moved to centre 
stage — and gives we feel far too central a role to thermodynamic issues.) 
In a similar vein, the energy equation 

V a (gV a )+p(V a V a )=0, (113) 

can be combined with the conservation equation 

V a (nV a )=0, (114) 

as follows: 

V a V a = — V a V a g = — p., (115) 

Q + p Q + p dr 

1 1 dr? 

V a V a = —V a V a n = — — , (116) 
n n dr 

where d/dr now refers to a material derivative along the flow. Eliminating the 
divergence we have 

J- ^ = i ^. (117) 
g + p dr n dr 

More formally, by invoking the barotropic condition we see that for any "fluid ele- 
ment" (i.e., tiny lump of fluid) we have 

^ = ^, (118) 
g+p n 

which can be rearranged as 

d{g/n) = -pd(l/n). (119) 

In terms of V, the "specific volume" of a fluid element, we have V oc 1/n and so 

d(g V) = -p dV, (120) 

which connects back to basic thermodynamics and again clearly verifies that in a 
relativistic setting g is the energy density, while n is the number density. 
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Finally, in view of these comments we can now write 



~- 2 V a V b + h° 



and somewhat simplify the acoustic metric to read 



;i2i) 



^^r'Tf^' 1 ^ <i22) 

/ \ 2/(d-l) (2 -\ 

Gah= {W^J \-^[Vo}a[Vo} b + h ab \. (123) 
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